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The control of shock noise or screech from a jet near a flexible structure is discussed. The pressure from the 
supersonic jet consists of a shock with spiral and flapping nonaxisymmetric modes superimposed on broadband 
response. This shock induces a nonlinear-nonstationary loading problem associated with acoustic wave generation 
and propagation coupled with structural vibration. Control of the shock is achieved by placing a ring at the nozzle 
lip oscillating at the shock fundamental frequency. The ring prevents the shock characteristics originating in the 
column of the shear layer from sustaining connection with the out-of-phase surface vibration. Shock-free flow is 
maintained over a large pressure ratio. The peak power spectral density of the pressure due to shock fundamental 
and harmonics is converted to broadband level. The experiment is motivated by consideration of aircraft structural 
fatigue and interior noise. 


L Introduction 

T HE noise generated by the exhaust of high-performance jet 
engines such as those used in the F-15 or B- 1 aircraft is a con- 
tinual concern. As a result, a method for controlling the shock noise 
or screech has been investigated. We demonstrate that a shock-free 
jet can be achieved by activating the nozzle lip in a feedback loop 
that interferes with the shock formation. This idea is used to suppress 
shock loading and reduce interior noise on a nearby aircraft fuse- 
lage structure in a laboratory. 1 Control of the shock at the nozzle lip 
means either 1) the shock is weakened or is converted into an expan- 
sion wave, all or in part, or 2) the vibrating lip attachment prevents 
the characteristics originating in the column of the shear layer from 
sustaining connection with an out-of-phase surface vibration. Shock 
from the jet induces high loading by the fundamental and harmonics 
tones at different locations on the structure. Controlling the shock at 
the nozzle exit has a greater advantage than controlling the structure. 

In the experiment, a model jet is used with a full-scale frame- 
stringer structure. As a result, the structure responds at a higher 
frequency compared with a full-scale nozzle placed at an equiva- 
lent geometrical distance. The jet flow is supersonic with an under- 
expanded shock. Active control of the shock at the nozzle exit is 
achieved by replacing the rigid lip with one oscillating at the shock 
fundamental frequency. 1,2 

Time series of the experimental data are used to unravel the com- 
plicated dynamics of the shock motion superimposed on the broad- 
band response using methods described by Abarbanel, 3 Abarbanel 
et al., 4 and Maestrello. 5 The pressure field develops into spirals and 
flapping nonaxisymmetric modes. The dynamics are captured by 
arrays of pressure transducers and used in the development of the 
shock control device. The structural loading due to the shock fun- 
damental and harmonics superimposed on the broadband response 
exhibits highly complex spatial chaotic deformation patterns of the 
type described by Dowell, 6 Nayfeh, 7 and Dowell. 8 The loading is 
sensitive to the initial condition due to the rotating and counterro- 
tating spirals randomly shed from the nozzle lip. Shock shells form 
between the shear layer and the lip. The tone is produced by the im- 
perfectly expanded jet. The literature refers to this as a screech tone. 9 
Hydrodynamics and structural stability analysis are the recognized 
tools that provide the loading description. 

Early experiments on shock jet noise by Powell 9 first recognized 
that the instability of the jet column due to the feedback loop be- 
tween fluid flow and sound is the cause of a powerful acoustic tone. 
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He proposed a simple mechanism for screech noise; he assumed 
that the sound waves moving upstream are of sufficient strength to 
govern the stability of the boundary layer near the exit. Later works 
by Westley and Woolley, 10 Norum, 11 Seiner, 12 Ponton and Seiner, 13 
Tam, 14 Morris et al., 15 and others indicated that the self-excited os- 
cillation involved the transfer of energy from one wave to another 
and that the spinning of the shock is associated with large-scale in- 
stability and sound. Tam presented theoretical arguments to show 
that the acoustic waves that radiate from the nozzle lip are con- 
fined to a relatively narrow frequency band if they are generated 
by the interaction between the large-scale instability waves and the 
quasi -periodic shock-cell structure of the jet. 

Section II of this paper describes the instrumentation and active 
controller. The tools used to analyze the dynamics response are 
described in Sec. III. Section IV describes the results, with Sec. 
IV.A describing the near-field pressure and shock control and Sec. 
IV.B describing the response of the structure with and without shock 
control. 

II. Instrumentation and Active Controller 

The experiment is carried out inside an anechoic chamber. 
The structure consists of six full-size aluminum aircraft-type pan- 
els, 45.72 x 17.71 x 0.1016 cm, with frame and stringers that are 
mounted on a rigid baffle and placed in the vicinity of the jet exhaust 
(Fig. 1). The baffle is inclined at an angle of 12 deg with the jet axis 
such that the jet plum will not interact with the baffle. 1 The struc- 
ture represents a section of a full-scale aircraft fuselage sidewall. A 
convergent nozzle with diameter D of 5.08 cm is used at a constant 
pressure ratio pot p a of 3 and stagnation temperature of 35°C, with 
corresponding Reynolds number U D/v of 1.3 x 10 6 and Strouhal 
number fD/U of 0.28, where p 0 and p a are the stagnation and 
atmospheric pressure, and /, U, and v are the shock fundamental 
frequency, jet exit flow velocity, and kinematic viscosity, respec- 
tively. The acoustic pressure radiated from the jet is measured by 
arrays of unvented, phase-amplitude-matched pressure transducers 
that are capable of dc response and placed in the vicinity of the 
nozzle exit. The vibration response of the structure is measured by 
two strain gauges, one located at the midpoint and the other on the 
right-side edge of the bottom panel (Fig. 1). These locations are 
sensitive to capture the passage of the shock during the rotation and 
flapping cycles. 

A new idea for controlling the shock at the exhaust is achieved by 
replacing the rigid nozzle lip with an oscillating one. The standard 
round nozzle and two views of the standard nozzle with an oscillating 
ring placed at the lip of the standard nozzle are shown in Fig. 2. The 
ring is coaxial with the jet and oscillates with the same frequency 
as the fundamental frequency of the shock. The ring is activated 
by a feedback loop from the pressure transducer that is placed in 
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the vicinity of the jet exhaust and supplemented by dc bias and 
a phase filter. 1,16,17 The pressure transducers monitor the output 
waveform in the vicinity of the exhaust. The ring is activated by 
two electromagnetic actuators. The pressure transducer at (—a, 0) 
provides the output signal A(t) = q(-a,0, t). The accelerometer 
on the ring provides the input signal //(O =m(0, 0, / - b ), where b 
denotes the time delay or phase shift. The difference Q{t) between 
input and output is used as a control signal as follows: 

CO) = *</)[/,■(/) - AO)] = K0)Q0) 

where K(t) is an adjustable amplitude. Also we can control the time 
delay b(t) of /; (r) with respect to A(t ). The control is introduced into 
the system input as a negative feedback ( K > 0). When the control 
is accomplished, the output signal A(t) is close to //(f), and, hence, 
the controlling force C(t ) becomes small. In the control stage, at 
least two control level adjustments in K and in the phase lag b are 
made by the feedback system. 

III. Data Analysis 

The tools used to analyze the dynamics of the response of the 
jet near-field pressure and panel response are explained. 1 - 5 Sound 
from higher-speed jets and the resultant structural responses have 
complex behaviors that are more accurately measured by nonlinear 
system response than by a stochastic system. Thus, the time histo- 
ries of the near-field pressure from the jet and the panel strain are 
obtained. From the time histories, the power spectral density, the 
phase portrait, and the probability distribution are computed. For a 
nonstationary signal q(t, x) t such as the pressure fluctuation p(t , x) 



a) Standard round nozzle 



b) Top view of standard round nozzle with oscillating lip 



c) Side view of standard round nozzle with oscillating lip 
Fig. 2 Nozzle. 


or panel strain s(/, x), the instantaneous power spectrum at instant 
T is defined by 


P(L T) = 


1 

2n 


L 


T+I/2 


T -Ifl 


exp(i2n ft)q(t, x) df 


where T is chosen so that the experimental run contains the interval 
T - 7/2, T + 1/2 for a sufficiently large I. The probability density 
of q(T, x) is denoted by Q(r , T), where 


<? or 




T + l/2 


q(t,x) dt 


T-l/2 


Q(r, T) = — prob [q(T, x ) < r] 
dr 
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In chaotic dynamics, searching for a low-dimensional character- 
ization of the system is of great interest. Let q(t,x) be a mea- 
sured temporal signal or time series at position x that is em- 
bedded in a d-dimensional phase space by a time delay r. The 
set Z(/) = [zi Zd(01 is regarded as a trajectory in the 

d-dimensional phase space. The distance between two points Z(/,) 
and Z(tj) is given by d xh and for a small e > 0, let N d (n, e) be the 
number of pairs of points with distance d tj < b. Then the correlation 
sum C d (e) and the correlation dimension Z)(d), for given d, are 
defined by Grassberger and Procaccia 18 as 


C d (e) = lim 


2N d (n,e) 


oo n(n — 1) 


D d = lim 

E-* 0 


log C d {e) 
logs 


For computation, the parameters r and d must be chosen not too 
large, and the correlation dimension D d is estimated by 


D d & 


log C d (n, b) 
log £ 


for a sufficiently small s and large n. The estimated dimension D is 
taken as the asymptotic value of D d as the embedding dimension d 
increases. 

Given the estimated dimension D, the Lyapunov exponent, which 
is one of the most important characteristics in a dynamic system, 
can be computed approximately. Several methods exist for comput- 
ing the Lyapunov exponents. 1819 The Eckmann-Ruelle method 20 
is used herein. 

Consider Z(/) as the trajectory of a dynamic system in the phase 
space of dimension d=D\ obtain the tangent (linear) maps 7J = 
I } 2, . . . f k of this reconstructed dynamical system by a least-squares 
fit; decompose 7} into an orthogonal matrix ft and an upper triangu- 
lar matrix by T\ - Q\R\ and 7} ft _i = ft/?,, fori >2; and compute 
the Lyapunov exponents as 


\j = lim[l/(A: — 1)] log|(/?* - 1 , . . . , /? 2 /?i)/il 

for i = l,2,...,D. For details see Maestrello 5 as well as the 
Eckmann-Ruelle algorithm given by Conte and Dubois. 21 


IV. Results 

The time series of the pressure is used to investigate the dynamics 
of the convected and radiated pressure from the jet with an oscillating 
shock. The investigation is carried out in two parts: the first contains 
the near-field pressure with shock and shock control, and the second 
contains the corresponding uncontrolled and controlled responses 
of structure. 


A. Near-Field Pressure and Shock Control 

The measured real-time near-field pressure p(r)> the computed 
power spectral density P(f> T), and the phase plot of the computed 
p(t) vs pit) are shown in Figs. 3a-3c. The real-time pressure near 
the instant T is used to evaluate the instantaneous plots of spectrum 
and phase. The shock produced by the underexpanded jet consists 
of a fundamental and two harmonic tones or screech, /„ 2f s , and 
3 f s . The shock fundamental dominates the broadband spectrum by 
as much as 40 dB. Figure 4 shows the real-time pressure at 0 = 0, 
90, 180, and 270 deg in the same location as Fig. 3 over extended 
time. We see that the envelop of the pressure fluctuation is nearly pe- 
riodic; the amplitude envelope is modulated over time. The pressure 
describes the position of the shear layer at any instant and is bounded 
by an envelope formed by the amplitude of the curve. These rapid 
changes overshoot the curve of the steady state. The oscillation rises 
above the starting level; after reaching maximum, the amplitude de- 
creases to the starting value. Azimuthally, there is a good correlation 
between data at two diametrically opposite points. We can trace the 
occurrence of a peak (or valley) at different angular locations to 
illustrate the asymmetry of the shock wave from the jet exit. This 
occurrence suggests a flapping motion of the shear layer. The phase 
plot is skewed, indicative of rotation in the azimuthal plane. In the 



/, sec 

a) Near-field pressure 



b) Power spectral density 


pit) 


pit) 

c) Phase plot 

Fig. 3 Pressure about nozzle exit for standard nozzle with c^> = 90 and 
0 = 0 deg. 


0, deg 



t, sec 

Fig. 4 Pressure about nozzle exit for standard nozzle with (j> = 90 , 0 = 0, 
90, 180, and 270 deg. 
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a) Near-field pressure 



/, kHz 

b) Power spectral density 



pH) 

c) Phase plot 

Fig. 5 Pressure about nozzle exit for standard nozzle with oscillating 
lip with <(> = 90 and 0 = 0 deg. 

meridian plane, 9 = 0 deg, previous measurements show the winding 
and unwinding motion of the jet column going through amplitude 
and bandwidth changes described in Figs. 5 and 6 of Ref. 1 . These 
changes are due to the spreading of the shear layer as the veloc- 
ity varies downstream. The slow drift of the wave modes could be 
attributed to the convection rotation and counterrotation effects, a 
result of the helical modulation of the pressure field, discussed also 
by Ponton and Seiner 13 and Westley and Woolley. 10 Tam 34 deter- 
mined the screech frequency by using hydrodynamic stability using 
multiple scales. This result is in agreement with the experiment. 

Control of the shock noise is achieved by replacing a rigid lip with 
an oscillating one. The new lip consists of a ring coaxial with the 
jet and oscillating with the same frequency as the fundamental of 
the attached shock. The control system is similar to the one adapted 
to control the structure response due to shock impingement. 1,2 ' 22 
The real-time pressure near the nozzle lip is used to activate the 
feedback control loop (Figs. 3 and 4). The ring is activated by two 
electromagnetic actuators with amplitude and phase variation K{t) 
and b(t) so that the energy of the harmonics is shifted back into 
the fundamental of the shock. Then, the controller reduces the am- 
plitude of the fundamental mode of the shock by increasing the 
amplitude adjusting to the control signal. Control of the shock is not 
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f, sec 
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a) Real-time strain 
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b) Power spectral density 
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/, kHz 


10 




r 


d) Probability 

Fig. 6 Panel strain induced by sound of jet exhausting from standard 
nozzle. 

straightforward, and full control is not always achieved in the initial 
stage of control. Some foreknowledge of the temporal dynamics is 
required during the initial stage. Several tries are needed to achieve 
a systematic feedback control, that is, K(t) and b(t), applicable 
for all initial stages. Shown in Fig. 5 are the measured real-time 
pressure p(r), the computed power spectral density P(f t T ), and 
the phase plot p(t) vs p{t) about the nozzle exit when the shock 
fundamental tone and harmonic tones are fully controlled. Control 
was applied to the fundamental frequency through a phase- and 
amplitude-tuning procedure in two stages to control the structural 
loading. Once control is achieved, the shear layer remains steady 
from flapping; thus, control is maintained. When comparing Fig. 5 
with Fig. 3, one notices a remarkable change in near-field pressure. 
The shock amplitude fundamental and harmonics reduce by several 
orders of magnitude to the broadband level. Similar changes occur 
to real-time plot, phase, and probability. Therefore, one can control 
the shock dynamics from the nozzle lip, rather than controlling the 
structure loading due to shock impingement. 

B. Response of Structure Without and with Shock Control 

The strain response exhibits broadband behavior with sharp, dis- 
tinct, high-amplitude spikes from the shock impingement by the 
fundamental and harmonics. The measured real-time strain s(/), 
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0 .4 .8 


r, sec 

a) Real-time strain 



/.kHz 

b) Power spectral density 




d) Probability 

Fig. 8 Panel strain induced by sound of jet exhausting from standard 
nozzle with oscillating lip. 


the computed power spectral density £(/, 7), the phase plot of the 
computed i(0 vs 5 (f), and the probability Q(r f 7) are shown in 
Fig. 6. The shock-induced oscillation dominates the response, the 
flapping of the shear layer induces random modulation as indicated 
by the real-time strain, the spectrum indicates the presence of the 
fundamental and two harmonics, the phase indicates rotation of 
the shear layer, and the nonsymmetric probability comes from the 
nonsymmetry in the panel deflection. Thus, results strongly indi- 
cate that the response of the structure is the image of the pressure 
loading. 

The flapping of the shear layer maintains some regularity part of 
the time (Fig. 4); the effect on the panel is not regular as shown by 
the real-time plot in Fig. 6a. The response of the structure is further 
analyzed by using the time series of experimental data of the panel 
response to evaluate the correlation dimension D and Lyapunov 
exponent X by using the Grassberger and Procaccia algorithm. 18 
The dimension of the attractor is related to the number of degrees 
of freedom of the panel response, as described by Ruelle 19 and 
Screenivasan. 23 Figure 7 shows the computed exponents, X 2 , X 3 , 
and X 4 , vs the embedding time (D - l)r from the temporal strain of 
the panel response with shock loading (Fig. 6), where D is the corre- 
lation dimension and r is the time delay. The exponents are positive, 
indicating nonlinear response behaviors. Only the largest Lyapunov 
exponent is shown for the experimental data. The results give some 
information about the strain response of the structure that indicates 
temporal chaotic behaviors when forced by the sound from a super- 
sonic jet with shock. The findings satisfy the experimental objective 
in determining the response to be used in the control strategy. 

Once shock is controlled at the nozzle lip, the response of the 
structure remains temporally smooth over the broadband. Compar- 
ing the plot of Fig. 8 with that of Fig. 6, we see the remarkable 
changes in response. The response of the fundamental and harmon- 
ics due to shock impingement after control reduces to the equivalent 
broadband response level. Similar changes are observed in the time 
plot, phase, and probability. Therefore, controlling the shock at the 
nozzle lip rather than over the structure is of great advantage for 
several reasons. One method requires a single controller, whereas 
controlling the structure requires several controllers because the 
shock fundamental and harmonics radiate over a wide area and at 
different angles, 17,22 


V, Conclusions 

A technique is demonstrated for controlling the shock from an 
underexpanded supersonic jet. A modified nozzle with a vibrating 
ring becomes an efficient shock suppressor. The shock freejet is a 
quiet and more efficient jet. The response of the structure is also 
significantly reduced. The technique is not limited to a circular ge- 
ometry; other geometries can be used as weR. The controller can 
also be utilized when the nozzle pressure varies. 

Two kinds of instability are induced by the shock or screech: a 
rapid flapping or pulsation of the jet column to each side and heli- 
cal disturbance along the jet column with winding and unwinding 
motions. Activating the nozzle lip in a feedback loop prevents the 
shock wave from sustaining a connection with the lip surface; as a 
result, both the flapping and helical instabilities decay. 

The ability to control the shock is of much practical importance. 
The quality and consistency of the results suggest that a full-scale 
experiment would be relevant. 
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